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ABSTRACT 

Using astrometric VLBI observations, we have determined the parallax of the black hole X-ray 
binary V404 Cyg to be 0.418 ± 0.024 mas, corresponding to a distance of 2.39 ± 0.14 kpc, significantly 
lower than the previously accepted value. This model-independent estimate is the most accurate 
distance to a Galactic stellar-mass black hole measured to date. With this new distance, we confirm 
that the source was not super-Eddington during its 1989 outburst. The fitted distance and proper 
motion imply that the black hole in this system likely formed in a supernova, with the peculiar velocity 
being consistent with a recoil (Blaauw) kick. The size of the quiescent jets inferred to exist in this 
system is < 1.4 AU at 22 GHz. Astrometric observations of a larger sample of such systems would 
provide useful insights into the formation and properties of accreting stellar-mass black holes. 
Subject headings: astrometry — X-rays: binaries — stars: distances — stars: kinematics — radio 
continuum: stars — stars: individual (V404 Cyg) 



1. INTRODUCTION 

Stellar-mass black holes in accreting X-ray binary sys- 
tems provide unique probes of the physical properties 
and formation mechanism of black holes, as well as the 
physics of accretion and associated outflow. However, 
our understanding of black hole systems is hampered 
by the fact that distances to X-ray binaries are rela- 
tively poorly known, typically uncertain by 50% or more 
([jonker & Nelemans 2004). Standard spectroscopic es- 
timates are handicapped by the distortions induced in 
the secondary by its compact par tner, contamination b y 
emission from the accretion disc (jRevnolds et al.|[2008D . 
and uncertainties in estimates of interstellar absorption. 

Since these distance uncertainties are systematic rather 
than random, they cannot be reduced by making mor e 
observations. According to I Jonker fc Nelemans! (|2004D . 
these systematic uncertainties can artificially enlarge the 
difference between the quiescent X-ray luminosities of 
black hole and neutron star systems. If true this would 
undermine the claim e d evidence for an e vent horizon 
(iNaravan et all Il997t iMenou et all 119991 : iGarcia et aP 
1200 iD . Accurate distances are also required to inter- 



pret the measured proper motions of black hole X-ray 
binaries, by converting a measured proper motion into 
a physical speed, from which we can derive the pe 
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culia r velocity of the source (e.g ., iMirabel et al.l 
12002; 'Mhabcl & Rodrigud 120031 : iDhawan et all 
Miller- Jones et al. 2009). This can be used to place con- 
straints on any velocity kick the black hole might have 
received in its natal supernova, and hence on the for- 
mation mechanism of the black hole (jBrandt et al.1 119951 : 
iWillems et aT1l2005l: llVagos et al.ll20m 

The only direct, model-independent method of mea- 
suring distances is via trigonometric parallax, which has 
hitherto been impossible for black hole X-ray binary sys- 
tems, since they lie at distances of several kpc, and hence 
require sub-milliarcsecond astrometry to measure their 
parallaxes. Very Long Baseline Interferometry (VLBI) at 
radio wavelengths is currently the only technique avail- 
able for such high-precision astrometric measurements. 
However, black hole X-ray binaries spend the majority 
of their time in a faint quiescent state, where they are of- 
ten not detected in the radio band at current sensitivities. 
They become bright enough to be detected during out- 
bursts, but such outbursts are typically not sufficiently 
long or frequent and do not sample the Earth's orbit well 
enough to make a parallax measurement feasible. Dur- 
ing outbursts, the radio emission is often resolved at the 
high angular resolutions required for precision astrome- 
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Fig. 1. — Position of V404 Cyg at each of our 7 epochs of 
observation. Dashed line shows the best-fit proper motion and 
parallax. Over time, the source moves south-west on the sky (from 
top left to bottom right of the figure). Each position is labeled 
with the MJD of the mid-point of the observation. 

try, making it difficult to determine the true location of 
the binary system from epoch to epoch. Lastly, many 
black hole X-ray binaries are located in or close to the 
Galactic Plane, such that the scatter-broadening along 
the line of sight makes high-precision astrometry impos- 
sible in the centimeter waveband where current VLBI 
arrays have the highest sensitivity. 

Here we present High Sensitivity Array (HSA) obser- 
vations of V404 Cyg, the most luminous known black 
hole X-ray binary in qu iescence. With a mass func- 
tion of 6.08 ± O.O6M0 jC asares fc Charles! [1991 ■ the 
compact object is a dynamically-confirmed black hole, 
accreting m atter from a KO subgiant companion star 
(jCasares e t al. 1993). The system has a quiescent radio 
flux density of 0.3 mJy, with a flat spectrum in dicative of 
a self-absorbed compact jet ()Gallo et al.ll2005l ). although 
the jet s are not resolved at the angular resolution of the 
HSA ijMiller- Jones et all [2001 . The persistent, unre- 
solved radio emission makes this source a good target 
for astrometric observations to measure its parallax. 

2. OBSERVATIONS AND RESULTS 

We observed V404 Cyg every 3 months over the course 
of one year, using the HSA. Two of the observations 
were made at the times of maximum parallactic displace- 
ment in Right Ascension (R.A.), since the R.A. signal is 
larger, and, owing to the greater size of the array in the 
east-west dimension, has smaller error bars than that in 
Declination. We observed at a frequency of 8.4 GHz, in 
dual circular polarization, with an observing bandwidth 
of 64 MHz per polarization. To these four HSA observa- 
tions we added three arc hival datasets, two of which hav e 
already been reported ([Miller- Jones et all [20081 I2009D . 
and one in which the source was previou sly not signifi- 
cantly detected ( Mioduszewski et al.ll2008f l. but knowing 
the astrometric parameters of the source, we were able 
to detect it at the 5a level. A fourth archival epoch did 
not detect the source to a 3f7 limit of 0.36mJybeam~^. 
A summary of the observations is listed in Table [1] and 
shows the extent of the source variability. In only one 
of the epochs was the source brighter than the 3tT up- 
per limit of this archival dataset (in which none of the 
large dishes of the HSA were available), making the non- 
detection consistent with the expected source flux den- 
sity. 

In all cases, the observations were phase-referenced to 
a bright, nearby calibrator from the International Ce- 



lestial Reference Frame (ICRF) source list, J2025-f3343 
(|Ma et all Il998l ). located only 16.6' from the target 
source. This implies that our systematic errors, which 
scale with distance from the phase reference source, 
should be relatively small 30/ias; see below). We 
observed in 3-min cycles, spending 1 min on the phase 
reference source and 2 min on the target in each cycle. 
To maximize astrometric accuracy, all data below 23° 
elevation were discarded. We also discarded data taken 
during the short ti mescale flares seen in a n umber of the 
observations fe.g.. iMiller- Jones et al.|[2008i) . in case the 
assumed variations in the jet power responsible for the 
flaring events translated into positional offsets. To en- 
able us to estimate the systematic errors affecting the 
astrometry, we substituted every seventh scan on the 
target source with an observation of a check source, the 
ICRF calibrator J2023-h3153, located 1°.87 away from 
the phase reference source. Data were r educed a ccord- 
ing to standard procedures within AIPS ([Greise 3[2003), 
and the source position at each epoch was determined 
by fitting an elliptical Gaussian to the source. In no 
case did the source appear to be resolved. Since the 
first epoch of archival data assumed a different position 
for the phase reference source, we corrected all epochs 
to a common calibrator position of 20''25™10". 8421056, 
33°43'00". 2144316 (J 2000) before fitting for the proper 
motion and parallax of V404 Cyg. The measured posi- 
tions are shown in Figure (TJ together with the best-fit 
p arallax and p r oper m otion. 

iPradel et al.l (|2006l ) used simulations to derive an ap- 
proximate formula for the systematic errors affecting var- 
ious VLBI arrays. For our calibrator-target separation 
and source declination, the estimated systematic errors 
(for a mean value of the wet zenith path delay) are 26 fias 
in R.A. and 30 ^as in Dec. The fitted positions of the 
check source show an rms scatter of 13 /las in R.A. and 
31/zas in Dec. Since the check source is seven times 
further away from the phase reference source than the 
target, this provides a rigorous upper limit on the sys- 
tematic errors in position. As seen from Table [1] our 
astrometric accuracy is limited by signal-to-noise rather 
than systematic erro rs. Neverthel e ss, th e systematic er- 
rors estimated from iPradel et al.l (|2006( ) were added to 
the statistical errors in quadrature before using the sin- 
gular value decomposition (SVD) method (as detailed 
in Loinard et al. 2007) to fit for a reference position and 
proper motion in R.A. and Dec, and for the source paral- 
lax. The best-fitting astrometric parameters, taking the 
mid-point of the observations, MJD 54322, as the refer- 
ence date, were: 

ao = 20^24^03". 821432 ± 0^.000002 (J 2000), 

So = 33°52'01".90134 ± 0".00005 (J 2000), 

Ha COS (5 = —5.04 ± 0.02 mas yr"^ , 

Mi ^ -7.64 ±0.03 mas yr"\ 

7r = 0.418 ± 0.024 mas. 

The fitted parallax signal, with the best-fitting proper 
motion removed, is shown in Figure [21 The reduced- 
value of the fit is 3.2, which, with 9 degrees of free- 
dom, implies a > 95%-confidence res ult. A bootstrap 
data-stripping analysis, as detailed by IChatteriee et al.l 



A parallax distance to V404 Cyg 



3 





I 













2000 2200 2400 2600 2800 3000 3200 3400 3600 



MJD-51544 

Fig. 2. — Parallax signature of V404 Cyg in right ascension 
and declination. The best-fit position and proper motion have 
been subtracted from the measured astrometric positions. The 
grey points show the displacement of the check source from its 
weighted mean position (dotted line), shifted by 0.7 mas for clar- 
ity, and scaled by a factor 1/7 to correct for its relative distance 
from the phase reference source. This illustrates the scale of the 
systematic errors. 

(|2009l ). confirmed the validity of these results, giving me- 
dian values differing by less than Itr from those found by 
the SVD fit, albeit with slightly larger error bars owing 
to the data stripping. Our fitted parallax corresponds to 
a source distance of 2.39 ± 0.14kpc. 

3. DISCUSSION 

3.1. Extinction-based distances 

Our fitted source distance is significantly closer 
than the best previous distance estimate of 4.0l^;2 ^pc 
ponker fc Nelemansll2004 ). However, that estimate was 
derived using an assumed interstellar extinction of Ay = 
3.3, and di d not take into ac count uncertainties in the 
extinction. iHvnes et al.l ()2009[ ) fitted the full multiwave- 
length spectral energy distribution with a template spec- 
trum for a KO IV star and found a formal best fit red- 
dening of Ay — 4.0. Indeed, we find that increasing Ay 
to close to 4.0 brings the estimate of Jonker & Nelemans 
([2004) down into agreement with our results. This high- 
lights the importance of using accurate extinction values 
when estimating source distances. Given the typical un- 
certainties in the extinction towards black hole soft X-ray 
transients (jjonkcr & Nelemans 2004), we conclude that 
this issue is likely to affect the majority of such systems, 
rendering their distance estimates also uncertain. 

3.2. Source luminosity 

Our measured distance is significantly closer than 
the 3. 5-4.0 kpc commonly assu m ed for the source 
(e.g., [Jm ler fc Nelemans' '2004"; 'Gal lo et al.l l2005l : 
iBradlev et al. 2007; Corbel et al. 2008), reducing its 
luminosity by a factor of 2.2-2.8. The quiescent 0.3- 
lOkeV unabsorbed fiux of 1.08 x 10^-^^ ergcm"^ s~^, 
derived using a powe r-law model for the X-ray spectrum 
(|Bradlev et al.l l2007f l. im plies a q uiesce nt luminosity 
of ~ 7 X lO^^ergs"^ iMakinol ()1989[ ) measured a 
maximum flux of 17 Crab with Ginga during the 1989 
outburst of V404 Cyg, implying a 1-70 keV luminosity 
of 7.0 X lO^^ergs"^, which is of order 0.5I/Edd for the 
black hole mass of 12 ± 2Mq derived for V404 Cyg 
(jShahbaz et al.lll994f ). Thus the system was not super- 
Eddington during the 1989 outburst. Furthermore, in 
reducing the source radio and X-ray luminosities, our 
new distance will reduce the scatter in the radio/X-ray 



correlation (| Gallo et al.l I2003D . bringing the points 
measured for V404 Cyg into better alignment with those 
of GX 339-4. 

3.3. Peculiar velocity 

With a more accurate distance and proper motion for 
the source, and the updated Galactic rotation parameters 
(i?f7i = 8.4 kpc, 6(71 = 254kms~M gi ven bv iReid et all 
(|2009| ) (but note iMcMillan fc Binnevl l2009l) we can re- 
visit the analysis of iMiller- Jones et al.l (120091 ) to derive a 
more accurate peculiar velocity of 39.9±5.5kms~^. This 
error bar now includes the uncertainties in both the dis- 
tance and in all three space velocity components. This 
new value is significantly lower than the best-fitting pe- 
culiar velocity of 64kms ~^ derived for a 4 kpc distance 
(iMiller- Jones et al1l2009( ). and can easily be achieved via 
a Blaauw kick, such that no asymmetric supernova kick 
is required. The component of the peculiar velocity in 
the Galactic Plane is 39.6kms~^, which, compared to 
the expected velocity dispersion in the Galactic Plane 
of 18.9kms"HMignard 2000) for the likely F0-F5 pro- 
genitor of the donor star, is a 2.1cr result, implying a 
probability of only 0.038 that the peculiar velocity is a 
res ult of Galactic velocity d ispersion (for more details, 
see IMiller- Jones et aTl l2009l) . We therefore find it most 
probable that the peculiar velocity arises from a natal su- 
pernova kick, the magnitude of which is consistent with 
recoil due to mass loss (a Blaauw kick), with any addi- 
tional asymmetric kick being small. 

3.4. Size constraints 

'Mil ler- Jones et all ()2009[ ) placed an upper limit of 
1.3 mas on the source size, corresponding to a physical 
size of < 3 AU at our new distance. While we did not re- 
solve the source, our highest-resolution data, taken with 
a global VLBI array at 22 GHz, constrain the source size 
to < 0.6 mas, a physical size of < 1.4 AU. We can there- 
fore place a lower limit of 10^ K on the brightness tem- 
perature of the source. From the high brig htness tem- 
perat ure, the observed flat radio spectrum (jCallo et al.l 
|2005| ) characteri stic of a steady, partially self-absorbed 
conical outflow (|Blandford fc Konigl' ll979l ). and the lo- 
cation of the sou rce on the radio/X-ray correlation of 
I Gallo et al.l (|2003l ). at the high-lumi nosity end of w hich 
jets have been directly resolved (|Stirling et al.ll2001[ ). we 
infer that the observed radio emission is likely to arise 
from compact, steady synchrotron-emitting jets. 

Our upper limit on the jet size can be compared to 
that of the resolved jets in Cygnus X-1 (jStirhng et al.l 
1200 It ), with an angular size of 15 mas at 8.4 GHz. Since 
the jet size scales inverse ly with observing frequency 
(|Blandford fc Konig f 19"7^. this implies 5 mas at 2 2 GHz 
for a physical scale of 10(d/2kpc) AU. Since Heini (|2006D 

8/17 

found the jet length should scale as Lj where L^, is the 
jet luminosity, then the difference in radio luminosities 
of the two sources suggests a size scale of ^ 1.5 AU for 
V404 Cyg at 22 GHz. Thus we are beginning to probe 
the angular scales on which the jets may be resolved. 

The scatter in the residuals after our parallax flt 
(O.lOmas in R.A. and 0.14mas in Dec.) can also help 
constrain the size and stability of the jets, since any 
motion of unresolved emitting knots along the jet will 
change the brightness distribution and hence the mea- 
sured centroid position. In particular, a comparison of 
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source positions measured during flares with those from 
the adjacent non-flaring data from which the proper mo- 
tion and parallax were derived will constrain the changes 
in jet morphology during flaring events. If the magnitude 
of the flares correlates with jet power and hence jet size, 
we might expect brighter flares to show the largest po- 
sitional shifts. For our brightest flare, in which the flux 
density increased by a factor of 3, the 3ct upper limit to 
the shift in the centroid position is 0.3 mas. This con- 
strains the emitting region during the flare to be at an 
angular separation of < 0.45 mas (a physical separation 
of < 1 AU) from the surface in the steady jet where the 
optical depth is unity. While such small positional shifts 
imply that minor flaring events should not destroy the 
parallax signal, if they can be better constrained dur- 
ing larger flares or with a higher-sensitivity array, such 
positional shifts could in future provide a method for de- 
termining the size scale of quiescent jets even in cases 
where they cannot be directly resolved. 

3.5. Wider implications 

These observations demonstrate that it is feasible to 
measure a parallax distance to a sufficiently bright qui- 
escent black hole X-ray binary. We have shown that the 
radio jets inferred to exist in the quiescent state of X- 
ray binaries do not hinder the astrometry sufficiently to 
render a parallax measurement impossible. The advent 
of the EVLA and the ongoing VLBA sensitivity upgrade 
will allow us to extend this method to a number of sys- 
tems with fainter quiescent jets, as well as opening up the 
possibility of detecting gy rosynchrotro n emission from 
low-mass donor stars (e.g.. lGudell2002f ). With accurate 
distances to a sample of quiescent black hole and neu- 
tron star systems, it will be possible to quantitatively 
assess the claim that black holes are less luminous than 
neutron star systems (iNaravan et al]|1997t iMenou et al.l 
119991 : iGarcia et al.l l2001h . Accurate outburst luminosi- 
ties for a sample of black holes will help quantify the fac- 
tor by which such systems can exceed their Eddington 
luminosities, with implications for the nature of Ultralu- 
minous X-ray sources (ULXs) in external galaxies. 

With a sample of black hole proper motions and par- 
allax distances, it will be possible to derive accurate pe- 
culiar velocities to probe the formation mechanisms of 
black holes of different masses. The three main chan- 
nels are direct collapse (giving no peculiar velocity), or 
via a supernova and fallback onto the proto-neutron star, 
with either a relatively small recoil kick for a symmet- 
ric supernova, or a potentially large kick from an asym- 
metric explosion. In the fallback case it is currently 
unclear whether and how the velocity scales with the 
amount of matter falling back on the neutron star. A 
sample of black hole peculiar velocities would provide 
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much-needed observational constraints, and help to de- 
termine whether there is a mass cutoff between black 
holes forming via direct collap se and those which form 
via a supernova and fallback (jFrver fc Kalogeral[200l[) . 
Finally, one of the methods currently employed to mea- 
sure black hole spins, via fitting of the X-ray spec- 
trum during the thermal-dominant stat e of the source 
during outburst (|McClintock et al]|2006[ ) requires accu- 
rate distances. Better distance constraints will allow 
for improved comparisons with spins derived via fitting 
the g eneral relativistic d istortions in disk reflection lines 
(e.g.. iMiller et alll2009r ). providing further insights into 
the generation and distribution of black hole spins. 

4. CONCLUSIONS 

We have used HSA observations to measure the most 
accurate distance to a black hole to date, via the method 
of trigonometric parallax. The distance of V404 Cyg is 
2.39±0.14kpc, signiflcantly lower than the previous best 
estimate of the source distance, and demonstrating that 
uncertainties in the interstellar extinction can introduce 
significant errors into distance estimates. With the new 
distance, we can determine that V404 Cyg did not exceed 
its Eddington luminosity during the 1989 outburst, and 
that it is likely to have formed in a natal supernova, with 
the peculiar velocity being consistent with recoil from the 
ejected mass, and any asymmetric kick being small. We 
have further constrained the size of the quiescent jets in 
the source to < 1.4 AU. These observations demonstrate 
the feasibility of measuring parallax distances to black 
hole X-ray binaries using VLBI techniques, which could, 
with a larger sample of systems, firm up evidence for 
black hole event horizons, and provide new insights into 
black hole spins, the mechanism of black hole formation, 
and the nature of ULXs. 
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TABLE 1 

ASTROMETRIC OBSERVATIONS OF V404 CyG 



Project 


Date 


MJD^ 


Array 


Frequency 


Bandwidth 


Flux density'' 


R-A.-^ 












(GHz) 


(MHz) 


(mjy beam^^) 


20'' 24™ 


33°52' 


BM117 


1999 Apr 10 


51278.60(25) 


VLBA 


15.353 


32 


< 0.36=^ 






BM231 


2005 Sep 16 


53629.17(21) 


VLBA, VLA 


15.365 


32 


0.33 ± 0.05 


03^822182(5) 


01".91599(9) 


BG168 


2007 Dec 02 


54436.84(09) 


VLBA, GBT, VLA 


8.421 


32 


0.28 ± 0.03 


03^.821266(5) 


01". 89861(22) 


GM064 


2008 Jun 01 


54618.42(25) 


VLBA, GBT, VLA, EVN« 


22.220 


64 


0.13 ± 0.03 


03". 821125(2) 


01". 89549(8) 


BM290 


2008 Nov 17 


54787.98(09) 


VLBA, GBT 


8.408 


64 


0.44 ± 0.03 


03^.820888(3) 


01". 89138(9) 


BM290 


2008 Feb 15 


54877.76(11) 


VLBA, GBT, VLA 


8.408 


64 


0.17 ±0.02 


03^.820826(5) 


01". 88959(12) 


BM290 


2009 Apr 26 


54947.57(11) 


VLBA, GBT, VLA 


8.408 


64 


0.33 ± 0.02 


03^.820777(3) 


01". 88813(8) 


BM290 


2009 Jul 03 


55015.38(10) 


VLBA, GBT, VLA 


8.408 


64 


0.13 ± 0.02 


03^820688(5) 


01". 88741(17) 



^ Specified Modified Julian Dates (MJDs) are for the mid-point of the observation, with the quoted uncertainties reflecting the observation duration 
^ Flux density measured after discarding short-timescale flaring events 

All positions are specified in J 2000 co-ordinates 
'' Source not significantly detected at this epoch; 3a upper limit given on flux density 

" Data taJcen with a global VLBI array, comprising the VLBA, GBT, phased VLA, seven dishes of the EVN (Cm, Eb, Jb2, Mc, Mh, Nt, Tr), and the Robledo antenna from the DSN 



